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Abstract. The grid of the models of radiative shock waves propagating through partially ionized hydrogen gas 
with temperature 3000K < T\ < 8000K and density 10~ 12 gmcm~ 3 < pi < 10~ 9 gmcm -3 is computed for shock 
velocities 20 kms -1 < U\ < 90 kms" 1 . The fraction of the total energy of the shock wave irreversibly lost due to 
radiation flux ranges from 0.3 to 0.8 for 20 kms" 1 < Ux < 70 kms" 1 . The postshock gas is compressed mostly due 
to radiative cooling in the hydrogen recombination zone and final compression ratios are within 1 < pN /pi < 10 2 , 
depending mostly on the shock velocity U\. The preshock gas temperature affects the shock wave structure due 
to the equilibrium ionization of the unperturbed hydrogen gas, since the rates of postshock relaxation processes 
are very sensitive to the number density of hydrogen ions ahead the discontinuous jump. Both the increase of 
the preshock gas temperature and the decrease of the preshock gas density lead to lower postshock compression 
ratios. The width of the shock wave decreases with increasing upstream velocity while the postshock gas is still 
partially ionized and increases as soon as the hydrogen is fully ionized. All shock wave models exhibit stronger 
upstream radiation flux emerging from the preshock outer boundary in comparison with downstream radiation 
flux emerging in the opposite direction from the postshock outer boundary. The difference between these fluxes 
depends on the shock velocity and ranges from 1% to 16% for 20 kms -1 < Ui < 60 kms -1 . The monochromatic 
radiation flux transported in hydrogen lines significantly exceeds the flux of the background continuum and all 
shock wave models demonstrate the hydrogen lines in emission. 
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1. Introduction 

Main indicators of shock waves in atmospheres of radially 
pulsating stars are enhanced hydrogen emission, double 
profiles of absorption lines and discontinuities in radial 
velocity curves. Most prominently these features are ob- 
served in W Vir QLebre fe Gillet 1992|), RV Tau ( pillct 
et al. 1989| ) and Mira type ( |Alvarez et al. 200fj| ) pulsat- 



(Hoefner & Dorfi 1997) it is nevertheless undoubted that 
shock-driven gas flows are responsible for formation of 
circumstellar shells. It must be noted that the postshock 
velocity near the photosphere should not necessarily be 
greater than the local escape velocity, because periodic 
passage of shock waves distends the stellar atmosphere 
and, under some conditions, may lead to the gradual ap- 



rauo oi trie aimospuere 
to the'sLellar radius in Lhese sLars is H/R > 10 _? and 
the linear theory of stellar pulsation admits the existence 
of progress ive waves propagating th rough the stellar at- 
mosphere (Unno 1965; Pijpers 1992). The running waves 
must inevitably transform into shock waves due to non- 
linear effects and this conclusion is corroborated by hy- 
drodynamic calculations of radial pulsations of low-mass 



proach of the gas flow to the escape velocity (Willson & 



in § varies. The laLio of Lhc aUiiosphcic scale height ffiU m ^ Furthermore; the shock amp ii tude tends to in- 



late-type gi a nts ([Wood 1974|; [Fadeyev fe Tutukov 1981 



Bowen 1988| ; [Fadeyev fe Muthsam 1990[ ). 

Though the role of shock waves as the principal mech- 
anism of mass loss from late-type stars is still disputable 
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crease while the shock wave propagates in a medium with 
decaying gas density. 

The most important problem concerning the shock- 
driven mass loss is the rate of the shock decay during 
its passage through the stellar atmosphere. Various sim- 
plifying assumptions lead to estimates of mass loss rates 
t hat differ fr om one another by many orders of magnitude 
( Wood 1979 ). Therefore, an analysis must be based on the 
self-consistent solution of the radiation transfer, fluid dy- 
namics and rate equations. Unfortunately, the description 
of the shock wave passage in terms of the Lagrangean 
approach cannot be applied, since it does not allow the 
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postshock relaxation zones to be explicitly resolved with 
respect to the spatial coordinate (Klein ct al. 1976, 197a ). 



Another approach is based on the assumption of steady- 
state gas flow. Applicability of this assumption for shock 
waves in stellar atmospheres is justified by the small width 
(including both the radiative precursor and the postshock 
relaxation zone) of the shock wave in comparison with the 
atmosphere scale height. 

Another necessary condition for applicability of the 
steady-state assumption implies that the characteristic 
time scale of the shock wave decay is idee <C d/U, where 
d and U are the width and velocity of the shock wave. 
In atmospheres of radially pulsating stars this condition 
is obviously fulfilled because emission lines indicating the 
shock wave ar e observed during the ma jor part of the pul- 
sation period ( Abt 1954 ; Willson 197€) and in so me stars 
this interval is as long as a half ( Fox et al. 1984 ) or even 
0.8 ( Gillet et al. 1985] ) of the pulsation period. 

Even in the framework of the steady-state assump- 
tion the problem of the shock wave structure in the par- 
tially ionized gas remains tremendously difficult, due to 
the strong coupling between the radiative precursor and 
the shock wake. Indeed, the occupation number densities 
of bound and free levels in both atoms and molecules 
strongly depend on the radiation field produced in the 
shock wake and at the same time must be determined as 
a solution of rate equations that are stiff on the relevant 
time scales. Fortunately, in stars with effective tempera- 
tures of T c g > 3000K, the gas is mostly in the atomic 
state and the principal dissipative processes in the post- 
shock region are excitation and ionization of atoms. 

In Paper I (Fadeyev & Gillet 1998) we proposed 



the method of global iterations that allows us to ob- 
tain the stable self-consistent solution of the equations of 
fluid dynamics, radiation transfer and rate equations for 
the steady-state plane-parallel shock wave pro pagating 



throug h the a tomic hydrogen gas. In Paper II ( Fadeyev 



& Gill :t 200C ) this method was employed for computa- 
tions of the structure of shocks with upstream velocities 
of 15 km s _1 < U\ < 70 km s^ 1 and for the temper- 
ature and the density of the unperturbed hydrogen gas 
of T\ = 6000K and p\ = 10~ 10 gm cm -3 , respectively. 
Results of these calculations are obviously insufficient for 
astrophysical applications since the temperature and the 
density in atmospheres of pulsating stars vary in wide 
ranges. 

Below we discuss the structure of shock waves propa- 
gating through partially ionized hydrogen gas with tem- 
perature and density ranges within 3000K < T\ < 8000K 
and 10~ 12 gm cm -3 < pi < 10~ 9 gm cm -3 , respec- 
tively. Basic equations and the method of their solution 
are described in Paper II. In comparison with our previ- 
ous work, we improved the treatment of radiation transfer 
and extended the total spectral range to 13 < log^ < 16. 
Furthermore, we employed the much more correct treat- 
ment of the spectral line radiation transfer which is of 
great importance for atomic level populations obtained as 
solution of rate equations. All calculations were done for 



the hydrogen atom with L = 4 bound levels and a contin- 
uum. 

2. General description of the shock wave structure 

In the comoving fluid frame the radiation-modified 



Rankine-Hugoniot relati ons are written ( Marshak 1958 
Mihalas fc Mihalas 1984]) as 



pU = Gq — m, 



mU + P K + P R = e 1 , 



1 * 2 

—mil + mh - 



U (E R + P R 



(1) 



(2) 



(3) 



where p is the gas density, U is the gas flow velocity, h is 
the specific enthalpy, P g is the gas pressure, Pr, Pr. and 
Pr are the radiation energy density, radiation flux and 
radiation pressure, respectively. 

In planar geometry the total energy flux C2 given by 
relation (||) is constant along the spatial coordinate X, 
so that while the parcel of gas passes the shock wave the 
decrease of the flux of kinetic energy Pk = \mU' 1 is bal- 
anced by changes in the enthalpy flux Ph = mh and the 
radiation flux Pr. Throughout the shock wave the flux 
U (Pr + Pr.) is significantly smaller than other terms of 
the left-hand side of relation (||) and can be ignored. 

As in Paper II the shock wave model is represented by 
the comoving plane-parallel finite slab. The space coor- 
dinate X = is set at the viscous adiabatic jump which 
is treated as an infinitesimally thin discontinuous jump 
where hydrodynamic variables undergo an abrupt change. 
It is assumed that the space coordinate X of the gas ele- 
ment increases while the gas flows through the shock wave, 
so that in the preshock region X < and in the post- 
shock region X > 0. The coordinates of the preshock and 
postshock outer boundaries are denoted as X\ and Xn, re- 
spectively. The upstream radiation flux emerging from the 
preshock outer boundary is negative, whereas the down- 
stream radiation flux emerging from the postshock outer 
boundary is positive, that is, Pri < and Prat > 0. 

On the upper panel of Fig. |l|, we show the fluxes Pk, 
Ph and Pr, as a function of space coordinate X for the 
shock wave model with p\ = 10~ 10 gmcm~ 3 , T\ = 6000K 
and Ui = 60 km s _1 . For the sake of convenience we use 
the logarithmic scale along the spatial coordinate X, the 
preshock and postshock regions being represented by the 
left plot and the right plot, respectively. Below the upper 
panel of Fig. [I], we show as a function of X the fractional 
number density of hydrogen atoms in 2-nd state rc^/fin, 
the hydrogen ionization degree xu — n c /na, the electron 
temperature T e and the temperature of heavy particles 
(i.e. of neutral hydrogen atoms and hydrogen ions) T a , and 
on the lowest panel we give the plot of the compression 
ratio p/pi- 

We assume that at the preshock outer boundary A 1; 
the gas temperature and the gas density are the same as 
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Fig. 1. The structure of the shock wave with p\ = 
1(T 10 gmcm" 3 , Ti = 6000K and U\ = 60 krns" 1 . On the 
upper panel the solid lines, dashed lines and dot-dashed 
lines represent the flux of the kinetic energy Fk , the flux of 
the enthalpy Fh and the radiation flux Fr. On lower pan- 
els are given the fractional number density of hydrogen 
atoms in the 2-nd state r^/riH, the hydrogen ionization 
degree xn, the electron temperature T c (solid lines) and 
the temperature of heavy particles T a (dashed line) , and 
the compression ratio p/pi- 

those of the unpertubed medium that are denoted as T\ 
and pi, respectively. At the same time, the upstream ra- 
diation flux emerging from the shock wake affects both 
the bound level number densities and the ionization de- 
gree of hydrogen atoms. Therefore, the number densities 
of hydrogen atoms in i-th state rii and the number den- 
sity of free electrons n e at the preshock outer boundary 
are determined from the solution of the equations of sta- 
tistical equilibrium. In our calculations we tried to set the 



preshock outer boundary Xi as far as possible from the 
discontinuous jump in order to diminish the perturbing in- 
fluence of the shock wave radiation field on the preshock 
gas. 

The most prominent changes that the gas flow under- 
goes in the preshock region are due to the absorption of 
the Lyman continuum radiation. In Fig. |l| the radiative 
precursor is revealed by the growth of the enthalpy flux 
Fh and due to the corresponding changes of the radia- 
tion flux Fr. The length of the radiative precursor 8X p 
is proportional to the mean free path of Lyman contin- 
uum photons and increases with decreasing density from 
5X p ~ 10 3 cm at p Y = 1CT 9 gmcm" 3 to 8X p ~ 10 7 cm at 
pi = 10~ 12 gmcm -3 . 

Within the radiative precursor the divergence of the 
Lyman continuum flux is negative, that is, 



V- F Lyc = An J {q v - n u J v ) du < 0, (4) 
in 

where r\ v and K v are the monochromatic emission and ab- 
sorption coefficients, J„ is the mean intensity and vi is 
the threshold frequency for ionization from the ground 
state of the hydrogen atom. Thus, the gas temperature 
increases due to absorption of the Lyman continuum ra- 
diation while the gas parcel approaches the discontinuous 
jump. As a result, the gas density p just ahead of the 
discontinuity exceeds its unperturbed value pi, whereas 
the gas flow velocity U becomes smaller than U\. The 
change of the gas density p and the gas flow velocity U in 
the radiative precursor depends not only on the upstream 
velocity Ui but also on the temperature and density of 
the unperturbed preshock gas. At pi = 10~ 10 gm cm" 3 
and Ui — 75 krns -1 , the relative growth of the density 
within the radiative precursor ranges from Ap/pi « 0.03 
at Ti = 3000K to Ap/ Pl « 0.07 at Ti = 7000K. 

It should be noted that when the hydrogen gas is not 
completely ionized in the preshock region, the increase in 
the enthalpy flux is mostly due to ionization of hydrogen 
atoms. The rates of collisional transitions, both bound- 
bound and bound-free, are negligible in comparison with 
those of radiative transitions, so that the increase in Fh 
ahead of the discontinuous jump results mostly from the 
photoionization of hydrogen atoms. 

Across the discontinuity, the flux of the kinetic energy 
Fk undergoes an abrupt decrease which is exactly bal- 
anced by the same abrupt increase of the enthalpy flux 
Fh. Because population number densities per unit mass 
rii/ ' p and n c / p do not change across the discontinuity, the 
increase in the specific enthalpy is only due to the abrupt 
increase in the temperature of heavy particles T a and the 
electron temperature T c . Here we assume that the tem- 
perature of heavy particles just behind the discontinuous 
jump is given by the solution of the Rankine-Hugoniot 
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relations (@) - (g): 

™o. „ m \ 1 to£/~ r; 2 — 1 

T+ = T--Xn(T+-T-) + ~- 



5 n H /c rj 2 



9 F + — F~ 

z f_R f_R 

5 n^kU~ 



2(E+ + P+)r ] -i-(E R +p-) 



(5) 



whereas the electron temperature increases due to adia- 
batic compression: 



T+ = 



(6) 



where rj = p + / p~ is the compression ratio across the dis- 
continuous jump, 7 = 5/3 is the ratio of specific heats of 
the electron gas, and superscripts minus and plus refer to 
the variables defined ahead of and behind the discontinu- 
ous jump, respectively []. 

The mean free path of photons in the vicinity of the 
discontinuous jump is several orders of magnitude larger 
than that of hydrogen atoms and the radiation flux will 
remain constant to a high order of accuracy, even if we 
take into account the finite thickness of the viscous adia- 
batic jump. In our calculations, the space interval between 
two adjacent cell centers located ahead of and behind the 
discontinuity was set equal to 1 cm for all models. Thus, 
the last two terms on the r.h.s. of (||) can be omitted with- 
out loss of accuracy because they refer to the difference 
between the upstream and downstream values of the radi- 
ation energy density, radiation flux and radiation pressure. 

Excitation and ionization of hydrogen atoms behind 
the discontinuity are mostly due to radiative transitions. 
The timescale of photoexcitation of the lower bound lev- 
els is comparable with that of the electron temperature 
growth, whereas the photoionization of hydrogen atoms 
is much slower and begins with a substantial delay The 
increase in the hydrogen ionization leads to a small in- 
crease in the enthalpy flux and, correspondingly, to the 
small decrease in the flux of kinetic energy because the ra- 
diation flux remains almost constant. However, the most 
promiment process accompanying hydrogen ionization is 
the increase in the compression ratio p/ pi. For the model 
presented in Fig. [l], the compression ratio increases from 
pj p\ = 3.6 just behind the discontinuity, to pj p\ = 6.8 in 
the layers where the hydrogen ionization reaches its max- 
imum. 

At larger distances from the discontinuity the kinetic 
energy of the gas flow initially stored as the ionization en- 
ergy of the gas is converted into radiation. In Fig. [l] these 
layers are revealed as those of remarkable changes in 
and Fn, accompanied by the growth of the compression 
ratio pi p\. For models considered in the present study, the 
compression ratio reached at the postshock outer bound- 
ary ranges from 1 < p?j / pi < 10 2 depending mostly on 



1 In relation ([]) we corrected the typeset errors which ap- 
peared in relation (15) of Paper II. 



the upstream velocity U\ and less on the preshock tem- 
perature T\ and preshock density p\. 

The large increase in the postshock compression ratio 
rj = pj pi is due to the ionization of hydrogen gas and its 
subsequent radiative cooling. In order to roughly estimate 
the contribution of these effects we can write the energy 
conservation relation with the omitted radiation en- 
ergy density and radiation pressure terms as 



II 2 

- u \e 



E, 



inl 



Eh 



Ft, 



Rl 

pi m 



(7) 



P 



in 



where E t = |(tih + n c )kT/p is the specific energy in the 
translational degrees of freedom and E m is the specific 
energy of excitation and ionization of hydrogen atoms. 
The left-hand side of (Q) is written for the preshock outer 
boundary, whereas its right-hand side represents an arbi- 
trary layer of the postshock region. 

Relation (f7j) and the momentum conservation relation 
(0) can be rewritten as 



7? = 4 + 3- 



E h 



E t 



F] 



Rl 



rhEt 



(8) 



where we assumed that the specific energy in the transla- 
tional degrees of freedom of the preshock gas is negligible 
in comparison with that of the postshock compressed gas, 
that is, E t i/ri < E t . 

The first term in the right-hand side of (g) is the lim- 
iting compression ratio of the atomic hydrogen gas at the 
discontinuous jump. The second term in the right-hand 
side of (||) describes the compression of the postshock gas 
due to excitation and ionization of hydrogen atoms. For 
the model shown in Fig. |l|, the maximum hydrogen ion- 
ization degree is sh = 0.99 and the ratio of the specific 
energies is E- ln / E t ~ 2. Thus, the upper limit for the com- 
pression ratio beyond the layers of hydrogen ionization is 
r) re 10. 

However the most important contribution to the 
growth of the compression ratio comes from the third 
term in the right-hand side of (f7|). Indeed, for the model 
shown in Fig. |l| the ratio of the sum of radiative fluxes 
(note that Fjh is negative) to the flux of the transla- 
tional energy of the gas at the postshock outer boundary is 
(Fun — Fm) / (rhEpf) ss 18. Therefore, the maximum com- 
pression ratio at the postshock outer boundary is rj « 64. 

3. Effects of the preshock temperature 

Just behind the discontinuity, the temperature of elec- 
trons, T e , is substantially lower than the temperature of 
heavy particles, T a , and electrons acquire energy in elastic 
collisions with neutral hydrogen atoms and hydrogen ions. 
Because the elastic scattering crosssection of electrons and 
hydrogen ions is much larger than that of electrons and 
neutral hydrogen atoms, the efficiency of the energy ex- 
change between heavy particles and free electrons is very 



Yu.A. Fadeyev and D. Gillet: The structure of radiative shock waves. Ill 



sensitive to the hydrogen ionization degree, x^, ahead of 
the discontinuous jump. In particular, at x^ > 10~ 2 , the 
postshock equilibration between the electron temperature 
T e and the temperature of heavy particles T a is due to 
elastic scattering of electrons by hydrogen ions and the 
equilibration rate very rapidly increases with increasing 
Xjj. Also, in reverse, at ih < 10~ 2 , free electrons gain 
most of the energy from elastic collisions with neutral hy- 
drogen atoms and the rate of temperature equilibation is 
much slower. This process, however, is perceptible only in 
shock waves with weak ionization in the radiative precur- 
sor (that is, at upstream velocities U x < 30 kms -1 ), the 
preshock gas temperature being T x < 4000K. 

Thus, effects of the preshock gas temperature on the 
shock wave structure are mostly due to the preshock equi- 
librium hydrogen ionization. This is illustrated in Fig. ^| 
for shock wave models with T x = 4000K and T x = 8000K. 
In both cases, the preshock gas density and the upstream 
velocity are p x = 10 -10 gm cm -3 and U x — 40 kms -1 , 
respectively, so that the postshock equilibration of T c and 
T a is due to elastic scattering of electrons by hydrogen 
ions. As is seen from the upper panel of Fig. ||, the hydro- 
gen ionization degree just behind the discontinuous jump 
is j: h » 0.02 at T x = 4000K and x u w 0.46 at T x = 8000K. 
As a result, at a preshock gas temperature T x — 8000K, 
the rate of the energy gain by electrons in elastic collisions 
with hydrogen ions is larger by a factor of 60 than that at 
T x = 4000K. 

The increase in the preshock gas temperature T x is ac- 
companied by the stronger radiation flux emerging from 
the shock wave. This is, obviously, due to the higher hydro- 
gen ionization degree in the postshock region. The more 
gradual changes in the postshock electron temperature 
and of the ionization degree in shocks with higher preshock 
temperature lead to smaller compression ratios pj p\ (see 
the lower panel of Fig. ||). 

In Fig. |^ we give the plots of the final compression ra- 
tio pn I ' p x which is reached at the postshock outer bound- 
ary Xn in shock wave models with p x — 10~ 10 gmcm -3 . 
It should be noted that the final postshock gas density 
asymptotically tends to its limiting value with X —>■ oo 
and because in our study the shock wave model is rep- 
resented by the finite slab, the compression ratio pn / Pi 
gives only the lower estimate for the limiting value Poo/ p x . 
In our calculations we tried to set the postshock outer 
boundary as far as possible from the discontinuous jump. 
Unfortunately, the convergence of global iterations is very 
sensitive to the spatial coordinate of the postshock outer 
boundary and for too large values of Xn, the iterations 
diverge. Thus, the coordinate Xn was determined for each 
model from trial calculations as a compromise between the 
requirement of the postshock region and demands of con- 
vergence and accuracy. For this reason, the filled circles in 
Fig. ||, representing the models with the same value of T x , 
do not lie on smooth curves. 

At a fixed upstream velocity, U x , the Mach number 
M x = U x /a x decreases with increasing T x due to the tem- 
perature dependence of the adiabatic sound speed a x at 
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Fig. 2. The postshock hydrogen ionization degree xu, the 
rate of energy gain by electrons in elastic collisions with 
hydrogen ions per unit mass Q el , the electron temperature 
T e (solid lines) and the temperature of heavy particles T a 
(dashed lines); the radiation flux i 7 ^ and the compression 
ratio p/p x as a function of distance from the discontinuous 
jump in shock wave models with p x = 10~ 10 gm cm~ 3 , 
U x = 40 km s-\ T x = 4000K and T x = 8000K. 

the preshock outer boundary. In a hydrogen gas with neg- 
ligible ionization (p x = 10~ 10 gmcnr 3 , T x < 7000K) the 
adiabatic sound speed is a x oc \fTi and the compression 
ratio 

Plim/Pl = lM\ (9) 

corresponding to the isothermal shock wave is inversly 
proportional to T x . On the other hand, as is seen from 
Fig. ||, the final compression ratio Pn/pi at the postshock 
outer boundary also decreases with increasing preshock 
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Fig. 3. The compression ratio pn / pi at the postshock 
outer boundary as a function of upstream velocity Ui at 
pi = 1CT 10 gmcixr 3 and 3000K < Ti < 7000K. Each se- 
quence of models with the same preshock temperature T\ 
is represented by filled circles connected by dashed lines. 
Each curve is labeled with T\. The solid line shows the 
compression ratio of the isothermal shock wave propagat- 
ing through the gas with T\ = 3000K. 
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Fig. 4. The final compression ratio pn / pi at the post- 
shock outer boundary in units of the compression ratio of 
the isothermal shock wave pi; m //9i versus the preshock gas 
temperature T\ for models with p\ — 1CP 10 gmcm~ 3 and 
20 km s _1 < \]\ < 90 km s _1 . At fixed preshock temper- 
ature T\ the values of the Mach number M\ correspond 
to models with JJ\ = 20 km s _1 and JJ\ = 90 km s" 1 , 
respectively. 



temperature though this dependence is not as prominent 
as that of pn m / Pi- 

In Fig. |i| we show the ratio pn / pnm for shock wave 
models with p x = 10~ 10 gmcm" 3 , 3000K < T x < 7000K, 
20 km s _1 < U\ < 90 km s _1 . Each model on this plot 
is represented by a filled circle and for models with U\ = 
20 km s _1 and U\ — 90 km s _1 we give the values of the 
upstream Mach number Mi. 

At a fixed preshock gas temperature T\ the ratio 
Pn I 'piim decreases with increasing upstream velocity since 
the final compression ratio p^/pi grows with U± more 
slowly than piim/pi (compare dependencies of p^/pi and 
piim/pi shown for T\ = 3000K in Fig. ||). However 
the most interesting feature is that the increase in the 
preshock gas temperature T\ is accompanied by the grad- 
ual approach of the final compression ratio Pn/pi to its 
upper limit pn m /pi corresponding to the isothermal shock 
wave. This is due to the decrease of the Mach number Mi 
with increasing T x . 

4. Effects of the preshock density 

The two upper panels of Fig. |B| show the postshock hydro- 
gen ionization degree xu and the postshock temperatures 
T c and T a as a function of distance from the discontinuous 
jump for shock wave models with 10~ 12 gmcm -3 < pi < 
10~ 9 gmcm- 3 , Ti = 6000K and Ui = 50 kms" 1 . The 
two lower panels of Fig. || for the same models show the 
radiation flux Fr and the compression ratio p/pi- Note 
that we use the logarithmic scale for the plot of radiative 



flux and the deep minima of log |Fr| correspond to layers 
with Fr w 0. 

The main results of these calculations are as follows. 
First, the width of the postshock relaxation zone increases 
with decreasing preshock gas density p\. This is due to 
the fact that the mean collision time of particles is pro- 
portional to the gas density. Second, the radiative flux 
produced by the shock wave decreases with decreasing pi 
because the kinetic energy of the gas flow at the preshock 
outer boundary is proportional to p\. Third, more grad- 
ual ionization and recombination of hydrogen atoms at the 
lower preshock gas density pi leads to smaller postshock 
compression ratios. 

5. Radiative flux of the shock wave 

As seen from plots on the upper panel of Fig. [j], the 
upstream radiative flux Fm emerging from the preshock 
outer boundary and the downstream radiation flux Frjv 
emerging in the opposite direction from the postshock 
outer boundary are not equal in absolute value. This dif- 
ference increases with increasing upstream velocity. For 
example, at pi = 10~ 10 gm cm~ 3 and Ti = 6000K the 
ratio of these fluxes ranges from |Fri|/ Frn ~ 1-01 at 
Ui = 20 kms^ 1 to \F m \/F RN w 1.16 at Ui = 60 kms" 1 . 

However the much stronger asymmetry of the radiation 
field is revealed at frequencies v > V\ because the Lyman 
continuum flux is directed mostly upstream. The region of 
the effective transport of the Lyman continuum radiation 
encompasses the near vicinity of the discontinuous jump. 
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Fig. 5. The hydrogen ionization degree cch, the electron 
temperature T c and the temperature of heavy particles T a , 
the radiation flux Fr and the compression ratio p/pi as 
a function of distance from the discontinuous jump in the 
postshock region at T\ = 6000K and U\ = 50 kms -1 . In 
solid, dashed, dot-dashed and dotted lines are shown the 
dependencies corresponding to the preshock gas densities 
pi = 1CP 9 , 10~ 10 , 10~ n and 10~ 12 gmcm -3 , respectively. 
Each curve is labeled with log p\. In temperature plots, the 
lower and upper converging curves represent the electron 
temperature T c and the heavy particle temperature T a , 
respectively. 



Ahead of the discontinuity, this is the radiative precursor 
and behind the discontinuity the Lyman continuum flux 
becomes negligible where the total radiation flux changes 
its sign (i. e.: in the layers with Fr ss 0). The contribution 
of the Lyman continuum to the total radiation flux within 
this region rapidly increases with increasing upstream ve- 
locity. For example, for shock waves propagating in a gas 
with pi = 1CT 10 gmcur 3 and T x = 3000K, the ratio of the 
Lyman continuum flux to the total radiation flux increases 



0.2 





log |X (cm) 

Fig. 6. The ratio of the radiation flux in the Lyman con- 
tinuum to the total radiation flux F^ yc /Fn in the vicin- 
ity of the discontinuous jump of shock waves with an up- 
stream velocity U\ — 50, 60, 70 (solid lines) and 80 kms -1 
(dashed line). 



from F Lyc /F R « 10~ 2 at Ui = 40 kms" 1 to F Lyc /F R « 0.2 
at Ui = 60 kms -1 . The growth of the ratio Fl yc /Fr ceases 
at an upstream velocity U\ ~ 75 kms -1 , corresponding to 
almost full hydrogen ionization in the radiative precursor. 
At larger upstream velocities the fraction of the Lyman 
continuum radiation does not exceed « 40% of the total 
radiation flux, but the region of the effective transport 
of the Lyman continuum radiation becomes wider, both 
ahead of and behind the discontinuity. 

The plots of the ratio Fl yc /Fr for shock wave models 
with upstream velocities 50 km s _1 < Ui < 80 km s _1 
are shown in Fig. ||. It should be noted that the positive 
value of the ratio Fl yc /Fr implies that the total radia- 
tion flux and the Lyman continuum flux are both negative 
and are directed upstream. Negative values of this ratio 
imply that the Lyman continuum flux is directed down- 
stream, whereas the total radiation flux is still upstream. 
Oscillations of the ratio FLyc/FR in the postshock region 
are due to vanishing values of Fr. 

Thus, behind the discontinuous jump, the Lyman con- 
tinuum flux FLyc changes its sign at smaller distance than 
the radiative flux 



F(u < v-l) = / F v dv 



(10) 



transported at frequencies below the Lyman continuum 
edge v\. In Fig. ^ we depict the zones of the postshock 
region with upstream and downstream total radiative flux 
Fr in shock wave models with p\ = 10~ 10 gm cm~ 3 and 
Ti = 6000K. The zone with downstream Lyman contin- 
uum flux (Ft^c > 0) and upstream flux F(v < v\) is 
shown on this diagram by the shaded area. Space coordi- 
nates X[F(v < v\) =0] and A^(Fr = 0) coincide because 
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Fig. 7. The zones of upstream (Fr < 0) and downstream 
(Fr > 0) total radiative flux in the postshock region 
at pi = 1Q- W gm cm~ 3 and T x = 6000K. The shaded 
area represents the zone with downstream Lyman contin- 
uum flux and upstream flux F(y < v\). The dotted curve 
marked \ max sh gives the locus where the hydrogen ion- 
ization degree is a half of its maximum postshock value. 



most of the radiation is transported in these layers at fre- 
quencies below the Lyman continuum edge. 

In Fig. [7] we also show the space coordinates of the lay- 
ers where half of the ionized hydrogen atoms are recom- 
bined. The close values of X(Fr = 0) and X(^maxxji) 
clearly illustrate that the main mechanism of the shock 
wave energy dissipation is the ionization of hydrogen 
atoms. 

At upstream velocity U\ < 60 km s _1 , the post- 
shock gas remains partially ionized and the maximum 
ionization degree increases with increasing U\, while the 
space coordinate of the ionization maximum decreases. 
Correspondingly, the hydrogen recombination also occurs 
at a smaller distance from the discontinuous jump. 

At upstream velocities of U\ > 60 km s -1 , the hy- 
drogen is fully ionized and the width of the zone with 
xti ~ 1 increases with increasing U\. Correspondingly, the 
hydrogen recombines at larger distances from the discon- 
tinuous jump and the coordinate X(Fr — 0) increases 
with increasing U\. At full hydrogen ionization, the zone 
of the effective energy transport in the Lyman continuum 
spreads both upstream and downstream with increasing 
upstream velocity and at U\ > 80 kms -1 both the Lyman 
continuum flux and the total radiation flux change sign 
nearly in the same layers. 

Another interesting conclusion which follows from the 
diagram in Fig. is that the width of shock waves with 
partial postshock hydrogen ionization decreases with in- 
creasing upstream velocity U\, whereas the width of shock 
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Fig. 8. Monochromatic mean intensity J„ (solid lines) 
and the local Planck function B V (T C ) (dashed lines) as 
a function of frequency v for four distinct layers of the 
shock wave with p = 10~ 10 gmcm -3 , T\ = 6000K and 
U\ = 60 km s _1 . The upper pair of plots represent the 
preshock region, whereas other plots represent the post- 
shock region. 

waves with full postshock hydrogen ionization increases 
with increasing upstream velocity. 

6. Frequency dependent radiation field 

The radiation field produced by the shock wave is remark- 
ably non-equilibrium and any attempts to use an assump- 
tion of thermal equilibrium inevitably lead to large errors. 
In order to evaluate the degree of departure from thermal 
equilibrium, it is instructive to compare the monochro- 
matic mean intensity 3 V with the local Planck function 
Bv{T e ). In Fig. H the plots of J„ and B V (T C ) are shown 
as a function of frequency v for four distinct layers of the 
shock wave model with p = 10 -10 gm cm -3 , T\ = 6000K 
and U\ — 60 kms -1 . For the sake of convenience each pair 
of plots is shifted with respect to others. 

The upper pair of plots represents J v and B v (T e ) at 
X 10 4 cm. As is seen from Fig. [j], the photoionization 
of hydrogen atoms in these layers due to absorption of 
the Lyman continuum radiation becomes perceptible and 
this layer can roughly be considered as the boundary of 
the radiative precursor. Excess radiation in comparison 
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with Bv(T c ) at frequencies v > v\ is due to the fact that 
the optical depth of the radiative precursor decreases with 
increasing frequency v. 

The second and third pairs of plots represent the cell 
just behind the discontinuous jump (X = 0.5 cm) and 
the layer where the hydrogen ionization degree reaches its 
maximum (X w 40 cm). Between these layers, various 
relaxation processes redistribute the energy of heavy par- 
ticles among other degrees of freedom and the radiation 
field is furthest from equilibrium. 

The lower pair of plots represent the layers (X w 
1.2 • 10 5 cm) where most of the hydrogen atoms are re- 
combined and the gas temperature gradually approaches 
its unperturbed value T\ . The optical depth at frequencies 
of the Lyman continuum is so large that 3 V and B P {T C ) 
coincide at v > v\. 

All shock wave models demonstrate two common fea- 
tures. First, we see as excess of Lyman continuum radi- 
ation in the preshock region, due to the existence of the 
radiative precursor. Second, there is a lower mean inten- 
sity of J v in comparison with B v (T e ) 7 due to the small 
optical depth of the shock wave at frequencies v < v\ . 

The plots of the mean intensity J„ displayed in Fig. || 
reveal the presence of six spectral line features. These are 
Lya, Ly/3, Ly7, Ha, H/3 and Paa. In comparison with 
our previous work described in Paper II, we improved the 
treatment of the spectral line radiation transfer and con- 
sidered each line of frequency Vq within the frequency in- 
terval as wide as [vo{l — S),i>q(1 + 8)], where 8 = 5- 10~ 4 . 
The use of wide frequency intervals is necessary because 
of the strong Doppler broadening of line profiles behind 
the discontinuous jump. In particular, the radiation flux 
transported in spectral lines is not negligible and contri- 
bution of the spectral line radiation to the total radiation 
flux increases with increasing upstream velocity. The most 
important are the Lya and Ha lines. For example, in the 
shock wave model with p\ = 10 -10 gmcm -3 , T\ — 6000K 
and U\ = 60 km s -1 the upstream radiation flux in Lya 
and Ha lines emerging from the preshock outer boundary 
is w 2% of the total radiation flux. 

In Fig. ^| we show the plots of the radiation flux 
integrated within frequency intervals of the Lya, Ly/3, 
Ly7 (upper panel), Ha, H/3 and Paa (lower panel) lines. 
Comparing with plots of n^/ nn and ih displayed in 
Fig. U we can conclude that the spectral line radiation 
is produced mostly in the layers of hydrogen recombina- 
tion. Another interesting feature is that the upstream and 
downstream Lya radiation fluxes are nearly equal in ab- 
solute value, whereas the radiation flux in Balmer lines 
emerges mostly upstream. 

7. Radiative losses of the shock wave 

In order to estimate the irreversible loss occurring in the 
shock wave energy it is instructive to compare the val- 
ues of each term of the left-hand side of relation (||) at 
both boundaries of the slab. In Fig. ^ these quantities are 
shown for shock wave models with p\ = 10 -10 gmcm -3 , 
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Fig. 9. The radiation flux integrated over spectral line 
frequency intervals as a function of space coordinate X 
for the shock wave model with p\ = 10 -10 gmcm -3 , T\ = 
6000K and U x = 60 kms -1 . 

T x = 6000K and 20 km s -1 < U x < 80 km s -1 . As is 
seen by comparing the upper and lower panels represent- 
ing the preshock and postshock outer boundaries, respec- 
tively, we see that most of the kinetic energy of the gas 
flow is irreversibly lost, since the ratio of the fluxes of ki- 
netic energy at both boundaries of the slab ranges within 
10 2 < Fki/Fkn < 7 • 10 3 . At the same time, the increase 
of the enthalpy flux after the passage of the parcel of gas 
through the slab does not exceed 20%. Thus, most of the 
energy of the shock wave is lost due to radiation. 

The role of radiation in irreversible energy losses of the 
shock wave can be evaluated from comparison of the ra- 
diation flux emerging from the boundary of the slab with 
the total energy flux 62- In Fig. |ll| the plots of the ra- 
tio Frat/C2 at the postshock outer boundary are shown 
for two sequences of shock wave models with T\ = 3000K 
and T\ = 6000K. As is seen, the radiative losses increase 
very rapidly with increasing upstream velocity U\ , though 
there is also the dependence on the preshock gas temper- 
ature T\. 

8. Conclusion 

In this work we improved the treatment of radiation trans- 
fer and considered the structure of shock waves propa- 
gating through partially ionized hydrogen gas with den- 
sities from 10 -12 gmcm -3 to 10 -9 gmcm -3 . Our calcu- 
lations showed that while the unpertubed hydrogen gas 
is partially ionized, the preshock gas temperature affects 
the shock wave structure, mostly due to the preshock hy- 
drogen ionization because the postshock equilibration of 
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Fig. 10. The flux of kinetic energy Fk (solid lines), the 
flux of enthalpy Fh (dashed lines), the radiation flux Fr 
(dot-dashed lines), and the flux U(Er + Pr) (dotted lines) 
at the preshock (upper panel) and postshock (lower panel) 
outer boundaries of the slab as a function of upstream 
velocity U\ in the models with p\ = 1CP 10 gm cm~ 3 and 
Ti = 6000K. 

translational degrees of freedom is very sensitive to the 
number density of hydrogen ions. In general, a higher tem- 
perature in the preshock gas will lead to larger postshock 
ionization and to a stronger radiation field produced by 
the shock wave. At the same time, the increase in the 
preshock gas temperature is accompanied by the slight 
decrease of the postshock gas compression. 

The preshock gas density pi affects the shock wave 
structure mostly due to the fact that the preshock flux 
of kinetic energy as well as the radiation flux emerging 
from the shock wave are proportional to p\. Therefore, 
in the gas with lower density, the shock wave produces 
weaker radiation flux and less compression of the post- 
shock gas. Moreover, at lower gas density, the width of the 
shock wave increases. The width of the radiative precur- 
sor increases due to the smaller absorption in the Lyman 
continuum, whereas the width of the postshock relaxation 
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Fig. 11. The ratio of the downstream radiation flux to the 
total energy flux -Frat/C^ at the postshock outer bound- 
ary as a function of upstream velocity U\ for shock wave 
models with p\ = 10~ 10 gmcm~ 3 . Each curve is labeled 
with Ti. 



zone increases because of the proportionality of the rates 
of relaxation processes to the gas density. 

The increase in the shock wave velocity leads to higher 
rates of relaxation processes in the compressed gas behind 
the discontinuous jump and to higher postshock ioniza- 
tion. While the maximum ionization degree of the post- 
shock gas is less than unity, the width of the shock wave 
decreases with increasing shock wave velocity. However, 
when the postshock gas is fully ionized, the increase in 
the shock wave velocity leads to the extension of the post- 
shock zone of hydrogen ionization and, therefore, to an 
increase in the shock wave width. 

Strong radiative cooling of the postshock gas leads to 
compression ratios as high as p/pi ~ 10 2 . Such large in- 
creases in the postshock gas density obviously can favour 
the condensation of dust grains in outer atmospheres of 
radially pulsating late-type giants. 

The isothermal approximation undoubtedly overesti- 
mates the radiative losses of the shock wave energy by at 
least a factor of two for shocks with an upstream veloc- 
ity of Ui < 30 km s _1 and can provide good accuracy at 
upstream velocities of U\ > 80 km s _1 , that is at Mach 
numbers Mi > 10. 

Our calculations demonstrated that the monochro- 
matic radiation flux at frequencies of hydrogen lines sig- 
nificantly exceeds the flux of the background contin- 
uum, so that our models reproduce the emission spec- 
trum observed in radially pulsating late-type supergiants. 
Moreover, the fraction of the energy transported in hydro- 
gen lines is not negligible and can be as large as one or two 
percent of the total radiation flux because of the strong 
broadening of the line profiles behind the discontinuous 
jump. 
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In this paper, we did not discuss the details of the emis- 
sion line spectra because current results are insufficient for 
comparison with observations. The change in the velocity 
of the gas during its passage through the shock wave must 
inevitably influence the line profiles. Calculations of shock 
wave structure based on the solution of the transfer equa- 
tion which takes into account effects of the gas flow ve- 
locity gradient were beyond the scope of the present work 
and will be presented in the forthcoming paper. 
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